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C/3 

■ Abstract. We review recent important topics in quantized vortices and quantum 

turbulence in atomic Bose-Einstein condensates (BECs) . They have previously been 
studied for a long time in superfluid helium. Quantum turbulence is currently one 
Ch ' of the most important topics in low-temperature physics. Atomic BECs have two 

distinct advantages over liquid helium for investigating such topics: quantized vor- 
tices can be directly visualized and the interaction parameters can be controlled by 
the Feshbach resonance. A general introduction is followed by a description of the 
dynamics of quantized vortices, hydrodynamic instability, and quantum turbulence 
in atomic BECs. 



1.1 Introduction 

Bose-Einstein condensation is often considered to be a macroscopic quantum 
phenomenon because bosons occupy the same single-particle ground state 
below the critical temperature for Bose-Einstein condensation so that they 
have a macroscopic wave function (order parameter) ^{r,t) = |iZ'(r, i)|e*^'-'''*^ 
that extends over the entire system. Here, the absolute squared amplitude 
I'f'P — n gives the condensate density and the gradient of the phase 9{r,t) 
gives the superfluid velocity field Vs = {h/m)\79 with boson mass m as the 
potential flow. Since the macroscopic wave function should be single-valued 
for the space coordinate r, the circulation F = § Vg-dt for an arbitrary closed 
loop in the fluid will be quantized with the quantum k — hjm. K vortex with 
such quantized circulation is known as a quantized vortex. Any rotational 
motion of a superfluid is sustained only by quantized vortices. Hydrodynamics 
dominated by quantized vortices is called quantum hydrodynamics (QHD), and 
turbulence comprised of quantized vortices is known as quantum turbulence 
(QT). 

A quantized vortex is a stable topological defect that is a characteristic 
of a Bose-Einstein condensate (BEC). It differs from a vortex in a classical 
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viscous fluid in the following three ways. First, unlike a classical vortex that 
can have an arbitrary circulation, the circulation of a quantized vortex is 
quantized. Second, since a quantized vortex is a vortex of inviscid superflow 
it cannot decay by viscous diffusion of vorticity, which occurs in classical 
fluids. Third, the core of a quantized vortex is very thin, being of the order 
of the coherence length (i.e., only a few angstroms in superfluid ''He and 
submicrometer in atomic BECs). Since the vortex core is very thin and does 
not decay by diffusion, the position of a quantized vortex in the fluid can 
always be identifled. 

Since any rotational motion of a superfluid is sustained by quantized vor- 
tices, QT usually takes the form of a disordered tangle of quantized vortices. 
QT is currently the most important research topic in QHD, which is an area 
in the fleld of low-temperature physics. The turbulence in a classical fluid, 
which is known as classical turbulence (CT), has been extensively studied in 
a number of fields, but it is still not well understood [1]. This is mainly because 
turbulence is a complicated dynamical phenomenon that is highly nonlinear. 
Vortices may represent the key for understanding turbulence, but they arc 
not well defined for a classical viscous fluid. They are unstable and appear 
and disappear repeatedly. The circulation is not conserved and varies between 
vortices. Comparison of QT and CT reveals deflnite differences, which demon- 
strates the importance of studying QT. QT consists of a tangle of quantized 
vortices that have the same conserved circulation. Thus, QT can be easier to 
study than CT and it offers a much simpler model of turbulence than CT. 

Quantized vortices and QT have historically been studied in superfluid 
helium. However, the realization of Bose-Einstein condensation in trapped 
atomic gases in 1995 provided another important system for studying quan- 
tized vortices and QT. The existence of superfluidity has been confirmed by 
creating and observing quantized vortices in atomic BECs and a lot of ef- 
fort has been devoted to studying lots of fascinating problems. Atomic BECs 
have several advantages over superfiuid helium, the most important being that 
modern optical techniques can be used to directly control their properties and 
to visualize quantized vortices. 

In a weakly interacting Bose system at zero temperature, the macroscopic 
wave function ^{r,t) obeys the Gross-Pitaevskii (GP) equation [2,3]: 

ij^^^hll = (-|^V2 + U{r) + gMr, t)\' - M^{r, t). (1.1) 

Here, U{r) represents the external potential (trapping potential, obstacle po- 
tential, etc.), g = 'I'Kh^a/m denotes the strength of the interaction character- 
ized by the s-wave scattering length a, and fi is the chemical potential. The 
only characteristic length scale of the GP model is the coherence length. It 
is defined by ^ = ?i/ {^/2rng\\^'\) and gives the vortex core size. The GP equa- 
tion has been used to interpret many properties of BECs in a dilute atomic 
gas [4,5]. It can explain both vortex dynamics and vortex core phenomena, 
such as reconnection and nucleation. 
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This chapter is organized as follows. Section 1.2 describes the dynamics 
of quantized vortices. It is very difHcult to directly observe their dynamics in 
superfluid "^He and '^He, which is one important advantage of atomic BECs 
over superfluid helium. In Sec. 1.3, we discuss hydrodynamic instability and 
QT in atomic BECs. The last section presents the conclusions. 

1.2 Dynamics of quantized vortices 

Since the realization of atomic BECs, many researchers have formed vortices 
by stirring a condensate or phase engineering. The details of these studies 
have been reviewed in other articles [6,7]. This chapter describes some recent 
contributions to this topic. 

1.2.1 Dynamics of a single vortex and vortex dipoles 

To gain insight into diverse superfluid phenomena, it is essential to understand 
the dynamics of a single vortex or vortex dipoles (i.e., vortex-antivortex pairs) . 
A few important observations have recently been made. Neely et al. nucleated 
vortex dipoles in an oblate EEC by forcing superfluid flow around a repul- 
sive Gaussian obstacle generated by a focused blue-detuned laser beam [8]. 
The nucleated vortex dipole propagates in a BEC cloud for many seconds; 
its continuous trajectory was found to be consistent with a numerical sim- 
ulation of the GP model. Freilich et al. observed the real-time dynamics of 
vortices by repeatedly extracting, expanding, and imaging small fractions of 
the condensate to visualize the motion of the vortex cores [9] . They nucleated 
vortices via the Kibble-Zurek mechanism [10, 11] in which rapid quenching 
of a cold thermal gas through the BEC phase transition causes topological 
defects to nucleate [12]. This produces single- vortex precesses in the cloud, 
whose frequency is in good agreement with a simple theoretical analysis [13]. 
Freilich et al. also nucleated vortex dipoles and observed their real-time dy- 
namics. This new technique enables quantitative comparisons to be performed 
between experiment and theory for vortex dynamics [14, 15]. 

Vortex dipole generation has been investigated theoretically and numer- 
ically in some recent studies based on the GP model. Sasaki et al. studied 
vortex shedding from an obstacle potential moving in a uniform BEC [16]. 
The flow around the obstacle is laminar when the velocity of the potential is 
sufficiently low. When the velocity exceeds a critical velocity of the order of 
the sound velocity, the potential commences to emit vortex trains. The man- 
ner in which vortices are emitted depends on the velocity and the width of 
the potential. The first pattern is the emission of the alternately inclined vor- 
tex pairs in a V-shaped wake, as shown in Fig. 1.1(a). The second pattern is 
produced by sequential shedding of two vortices having the same circulation. 
The two vortices rotate about their center without varying their separation, 
forming a train of vortices similar to a Benard-von Karman vortex street, as 



4 



Makoto Tsubota^ and Kenichi Kasamatsu'^ 



shown in Fig. 1.1(b). For a wide potential with a high velocity, the period- 
icity disappears, as shown in Fig. 1.1(c). Aioi et al. subsequently proposed 



(a)? = 2.4,d/^ = 0.04 (b) T = 2.6, d / 4 = 0.05 (c) T = 3.0, d / 1 = 0.05 




Fig. 1.1. Density and phase distributions of a condensate past an obstacle potential. 
The velocity and potential width are (a) {v,d/£) = (2.4,0.04), (b) (2.6,0.05), and 
(c) (3.0,0.05), where v = v{lQ^m/ gnof^'^ and ^ = n{lQ^ /{mgnof^ with no the 
atom density without perturbation. The white arrows in (a) indicate the directions 
in which the vortex pairs move. The field of view is 6^ x 3^. [Sasaki, Suzuki and 
Saito: Phys. Rev. Lett. 104, 150404 (2010), reproduced with permission. Copyright 
2010 by the American Physical Society]. 

controlled generation and manipulation of vortex dipoles by using several 
Gaussian beams from a red (attractive potential) or blue (repulsive potential) 
detuned laser [17]. For example, when a red-detuned beam moves through 
a BEC cloud above a critical velocity, a vortex dipole nucleates at the head 
of the potential. In contrast, a blue-detuned beam generates vortex dipoles 
on both sides of the potential. Double beams can generate various kinds of 
dipole wakes depending on the velocity and the distance of the two beams. 
The trajectory of emitted vortices can be modified by another potential. 

1.2.2 Dynamics of vortices generated by an oscillating obstacle 
potential 

Recently, in the field of superfluid '*He and '^He, several groups have experi- 
mentally studied QT generated by oscillating structures such as wires, spheres, 
and grids [18]. Despite significant differences between the structures used, their 
responses with respect to the alternating drive have revealed some surprising 
common phenomena. The response was laminar at low driving rates, whereas 
it became turbulent at high driving rates. 

This strategy can also be applied to trapped atomic BECs. Fujimoto and 
Tsubota numerically investigated the two-dimensional dynamics of trapped 
BECs induced by an oscillating repulsive Gaussian potential [19,20] and found 
a strong dependence on the amplitude and frequency of the potential. Unlike a 
potential with constant velocity [16], an oscillating potential continually sheds 
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vortex pairs with alternating impulses, a typical example of which is shown in 
Fig. 1.2. The nucleated pairs form new vortex pairs through reconnections that 
move away from the obstacle potential and toward the condensate surface. The 
BEC cloud eventually becomes full with such vortices, as shown in Fig. 1.2(d). 
An oscillating potential is thus a useful tool for generating QT in a trapped 
BEC. 




Fig. 1.2. Reconnection of vortex pairs about an oscillating potential. The symbols 
- and + indicate vortices with clockwise and counterclockwise circulations, respec- 
tively. The black arrows indicate the direction in which the potential moves. The 
horizontal and vertical dimensions of (a), (b) and (c) are 27.0 and 14.0 micrometer, 
while those of (d) are 145 and 34.0 micrometer. When the potential moves to the 
right, it generates a vortex pair with an impulse in the same direction as the poten- 
tial. When the potential (a) (t—SlA ms) changes direction and moves to the left, (b) 
(f=38.3ms) it collides with the pair and (c) (f=43.2 ms) emits another pair. These 
two pairs reconnect to form new pairs that move away from the potential, (d) The 
BEC cloud eventually becomes filled with such vortices. [Fujimoto and Tsubota: 
Phys. Rev. A 83, 053609 (2011), reproduced with permission. Copyright 2011 by 
the American Physical Society]. 



1.2.3 Kelvin wave dynamics 

Kelvin waves are three-dimensional excitations along a vortex line. Kelvin 
waves play an important role in dissipation in superfluid helium at very low 
temperatures in which the normal fluid component is negligible [21]. Kelvin 
waves have been experimentally demonstrated in a trapped BEC [22] by ex- 
citing a collective quadrupole mode to a single- vortex state that subsequently 
decays to Kelvin modes by a nonlinear Beliaev process. This was supported 
by numerical analysis based on the Bogoliubov-de Gennes equation [23] and 
the CP equation [24] . Rooney et al. recently studied the dynamics of Kelvin 
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waves using the stochastic projected GP equation [25]. They showed that 
Kelvin waves can be suppressed by tightening the confinenient of the trap 
along the vortex line, which drastically reduces the vortex decay rate as the 
system becomes two-dimensional. This behavior is consistent with observa- 
tions of the decay of vortex dipoles [8]. 

1.3 Hydrodynamic instability and quantum turbulence 

Some theoretical ideas for achieving QT in trapped BECs have been proposed. 
The turbulent state may be generated during the formation of a condensate 
from a nonequilibrium non-condensed Bose gas by rapid quenching [26] . Con- 
versely, turbulence is expected to be generated during the destruction process 
from the equilibrium condensed state. This subsection describes recent studies 
on hydrodynamic instability and the resulting nonlinear dynamics that may 
lead to QT. These instabilities are formed by applying an external driving 
force to condensates or by complicated interactive phenomena between mul- 
ticomponent BECs (some of which are analogs of well-known phenomena in 
classical hydrodynamics) . 

1.3.1 Methods to produce turbulence in trapped BECs 

A major problem when investigating QT in atomic BECs is the difficulty in 

applying a dc velocity field in superfluid helium. Here, we summarize some 
proposals for and a recent experimental realization of QT generation in a 
trapped single-component BEC. 

Simple rotation 

An important way for generating vortices in trapped BECs is to rotate the 
external potential [27-29] . However, rotation alone cannot lead to QT because 
it generates an ordered vortex lattice along the rotational axis [30], which 
is the equilibrium state in the corresponding rotating frame. Nevertheless, 
Parker and Adams suggested the emergence and decay of turbulence in a BEC 
under a simple rotation, starting from a vortex-free equilibrium BEC [31]. A 
numerical simulation based on the energy-conserving GP equation suggests 
the existence of a turbulent regime that contains many vortices and high- 
energy-density fluctuations (sound field) on a route to the ordered vortex 
lattice. 

Two-axis rotation 

Since the above turbulence is generated during the ordering process, it is not 
steady turbulence. Kobayashi and Tsubota suggested performing rotations 
about two axes [32], as shown in Fig. 1.3(a). When the spinning and precessing 
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rotational axes are perpendicular, the two rotations do not commute and 
thus cannot be represented by their sum. This situation can be modeled by 
simulations based on the GP equation: 

where L is the angular momentum and the rotation vector is written as fl{t) = 
(fix, sin Qxt, f^z cos f2xt) with frequencies of fl^ and i?^ for the first and 
second rotations, respectively. The inclusion of these non-commuting rotations 
and phenomenological dissipation which is effective only at scales smaller than 
the healing length successfully generate steady turbulence [32]. 




Fig. 1.3. QT in atomic BECs. (a) A method for realizing steady turbulence in an 
atomic BEC subject to precession. A BEC is trapped in a weakly elliptical harmonic 
potential. A rotation is applied along the z-axis followed by a rotation along the x- 
axis. (b) Energy spectrum of steady QT obtained by two-axis rotation. The dots 
represent the numerically obtained spectrum for an incompressible kinetic energy, 
while the solid line i s the K olmogorov spectrum. Here, Rtf is the Thomas-Fermi 
radius and atrap = \fh[mjj is the characteristic length scale of the trap. 



Donnelly-Glaberson instability 

When both rotation and a linear velocity are applied to the system, the right- 
hand side of the GP equation Eq. (1.2) will have the term —V ■ p in the 
corresponding helically moving frame. For V = Vz the dispersion of the vortex 
waves (Kelvin waves) on a straight vortex line parallel to the z-axis behaves 
as w — Vkz and its frequency will become negative above a certain critical 
velocity. This is known as the Donnelly-Glaberson instability in superfluid 
helium and it amplifies Kelvin waves [33]. If the rotating BEC contains a 
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vortex lattice, the amplified Kelvin wave induces reconnections of adjoining 
vortex lines, eventually leading to a turbulent state [34]. 

Combined rotation and oscillating excitation 

As a method that enables better control of turbulence, Henn et al. introduced 
an external oscillatory potential to an *^^Rb BEC [35,36]. This oscillatory field 
induced a successive coherent mode excitation in a BEC. They observed that 
increasing the amplitude of the oscillating field and the excitation period in- 
creased the number of vortices and eventually lead to the turbulent state [36]. 
In the turbulent regime, they observed a rapid increase in the number of vor- 
tices followed by proliferation of vortex lines in all directions, where many 
vortices with no preferred orientation formed a vortex tangle. The oscillatory 
excitation (which mainly consists of oscillation, rotation, and deformation) 
nucleated vortices. However, it is still not known theoretically how the turbu- 
lence is generated. 



1.3.2 Signature of quantum turbulence 

It is important to determine whether the highly excited state is really QT. 
Several methods are used to identify the turbulent state; they are discussed 

in this subsection. 

In most numerical simulations, apart from observing a random configu- 
ration of vortices, the turbulent regime has been identified by checking if 
its incompressible kinetic energy spectrum obeys Kolmogorov's —5/3 law 
[1]. When the condensate wave function is written in the form if'(r,t) = 
i/n(r,t)e*^('''*\ the kinetic energy is expressed by the sum i^kin = Ek + Eg, 
where Ek ~ {h /2m) J dr\y^V9\'^ denotes the superfluid kinetic energy and 
Eg = {h^ /2m) J dr|Vv^|^ is the quantum pressure energy. The vector field 
^/nVO can be divided into incompressible (solenoidal) and compressible (irro- 
tational) components: ^/nVO = {,/nVey + {./nVe)", where V • {^/nVOy = 
and V X [y/nVOy = 0. Thus, the incompressible and compressible kinetic 
energies are defined by E]:" = {h^ /2m) J drli^VeY'"]"^ . They correspond to 
the kinetic energies in the vortices and the sound waves, respectively. Since 
the compressible and incompressible fields are mutually orthogonal, it follows 
that Ek = E'l+E^. The kinetic energy spectrum as a function of wave number 
k is defined by 



^ ^ 2m J 



sin 



(27r)3 



(1.3) 



such that E]^'^ = ejf{k)dk. The Kolmogorov law states that the incom- 
pressible energy spectrum obeys the power law e\{k) ~ k" where v = —5/3 
over the inertial range of k. For a trapped BEC, the inertial range that follows 
the Kolmogorov law is determined by the Thomas Fermi radius i?TF and the 
coherence length ^ = h/2mii [32] [see the right panel of Fig. 1.3]. 
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The structure of QT is reflected in the time dependence of the decay of the 
total vortex line density L after turning off the excitation that sustains the 
turbulence. Correlations of vortex tangles can be classified into two kinds [37] : 
correlated and uncorrelated tangles. In a correlated tangle, turbulent energy 
is concentrated in the "classical" length scale range that is larger than the 
intervortex distance I, where the correlated tangle exhibits a Kolmogorov 
spectrum. The energy is then transferred to smaller scales by a Richardson 
cascade and L decays as In an uncorrelated tangle, the turbulent energy 

is associated with a random vortex tangle with spacing I and it is concentrated 
in the "quantum" range of length scales smaller than I. Then, L decays sst~^. 
Thus, observing the scaling behavior of i on f provides useful information 
about QT. 

In addition. White et al. [38] showed that QT is characterized by a power- 
law behavior of the probability density function of the velocity field Vs{r) — 
{h/m.)V0{r), whereas classical turbulence obeys Gaussian velocity statistics. 
This non-Gaussian behavior originates from the singular nature of a quantized 
vorticity with a 1/r velocity field, which appears in the high- velocity region 
determined by the intervortex distance v ~ k/itI [39]. 

It is difficult to measure the above statistical or scaling properties ex- 
perimentally. Henn et al. observed another remarkable feature of turbulent 
condensates: suppression of aspect ratio inversion during free expansion af- 
ter turning the trapping potential off [36] . Despite the asymmetric expansion 
(from a cigar shape to a pancake shape) of a conventional BEC or isotropic 
expansion of a thermal cloud, the turbulent state exhibited a self-similar ex- 
pansion that preserved the initial aspect ratio. Although a quantitative theo- 
retical understanding of this effect has yet to be fully realized [40], it represents 
a remarkable new effect in the turbulent regime. 

1.3.3 Hydrodynamic instability in multicomponent BECs 

Multicomponent atomic BECs can be created in cold-atom systems with, 
for example, multiple hyperfine spin states or a mixture of different atomic 
species. Such systems yield a rich variety of superfiuid dynamics due to the 
intercomponent interaction. Two-component BECs are the simplest multi- 
component system. Schweikhard et al. experimentally investigated the vortex- 
lattice dynamics of two interacting and rotating condensates by transferring 
some of the initial population of ^^Rb BECs with a vortex lattice to its other 
hyperfine state via a coupling pulse [41]. They observed the ordering dynamics 
change from a triangular lattice structure to a stable square lattice through 
the transient turbulent regime. 

Recent experimental advances have provided more controllable ways to 
study the rich dynamics of two-component BECs. External potentials can be 
prepared that can act independently on both components; this enables initial 
conditions to be prepared that are suitable for studying a particular problem. 
In addition, the intra- and inter-component interactions can be tuned with 
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the help of the Feshbach resonance [42-44]. This allows phase separation to 
be performed and interface phenomena between two superfluids to be studied 
in a well-controlled manner. 

In the mean-field theory, a two-component BEC is described by macro- 
scopic wave functions 'Pi, where the subscript i refers to each component 
{i = 1,2). The Lagrangian for this system is given by 

L = J dr{Pi + P2- gi2\^i\^\^2\'') , (1.4) 

where 

= ^m*^ + ^w:v^w, - u,m' - (i.s) 

at Zrui 1 

with m,i and Ui being the atomic mass and the external potential of the ith 
component, respectively. The intra- and inter-component interaction param- 
eters have the form 

Qij = 2Ti'h^aij{ml^ + mj^), (1.6) 

where aij is the s-wave scattering length between the ith and jth components; 
we assume > in the following. For homogeneous condensates, the con- 
densates are miscible and immiscible when 511522 > 9i2 ^-iid 5ii522 < 5i2i ''c- 
spectively. Below, we review the characteristic hydrodynamic instability that 
occurs in each condition. 



Counter-superflow instability 

When two-component BECs coexist with a relative velocity, they exhibit dy- 
namic instability above a critical relative velocity [45]. This phenomenon is 
known as a counter-swper/ilow instability (CSI). Takeuchi et al. suggested that 
the nonlinear dynamics triggered by the CSI generates a binary QT composed 
of coreless vortices and that thus has a continuous velocity field [46,47]. Ham- 
ner et al. realized the CSI experimentally for quasi-lD geometry and observed 
the generation of shock waves and dark-bright solitons [48] . 

The CSI can be understood from the Bogoliubov spectrum for a system of 
a uniform two-component BEC with a relative velocity [45]. The functional 
derivative of J dtL, where L is the Lagrangian in Eq. (1.4), with respect to 

gives the GP equation. 

In this subsection, we assume f/j = and the miscible condition 5115(22 > 5i2- 
The wave functions !Pi = ^'jo in a stationary state can be written as 

iZ^io = V^e^^™*'^*"-"**'/'^ (1.8) 
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with the velocity Vi and the chemical potential /Xj of the ith component. 

Counter supcrflow occurs when Vi ^ V2- Wc consider a small excitation 
above the stationary state as = !Z'jo + ^^ii where we write the excitation of 
the wave functions S'Ps in the form 



5\I/ = g^ii™iVi-r-Hit)/h 



(1.9) 



By linearizing the GP equation (1.7) with respect to 6^i, we obtain the 
Bogoliubov-de Genne equations: 



2m 

21.2 



hk-Vi + guUio^ Ui - gaHioVi + gij^nmnjo {uj - vj) = Jioji^l,W) 

( "2m~ ~ ■ ^ + gunioj Vi - guriioUi - gij y/fiionjo {uj - vj) = -?ia4i-ll) 

Diagonalizing the eigenvalue equations (1.10) and (1-11), we obtain the Bogoli- 
ubov excitation spectrum. For simplicity, we set mi = m2 = m, gn = 322 = g, 
and nio = 7120 = no, and neglect the center-of-mass velocity of the two con- 
densates. The eigenvalue of Eqs. (1.10) and (1.11) then has the simple form: 

nW =e' + \h'qV^ ± Je^h'BV^ + AglX4, (1.12) 



where Vr = Vi — 1/2 is the relative velocity and we denote = eo(eo + 2.9no), 
eo = n^k'^/2m, and fc^ fcj^ + with k\\ = (fc • Vr)/Vr and k^ > 0. 
The system is dynamically unstable when the excitation frequency uj be- 
comes imaginary (i.e., cj^ < 0), which reduces to e_ < ?ifc||Vn/2 < e+ with 
e± = \J eo[eo + 2((? ± gvi)'nQ)\. This inequality determines the unstable region 
in wavenumber space (fcn, fcj,), where Vr exceeds the critical relative velocity 
Vc = V- with V± = 2^ygno/m^yi ± 312/3- The total momentum density 5 J 
carried by the excitation can be defined as 6 J = 5J\ + 5J2 with a change 
in the momentum density 5Jj = h,k{\ui\^ — \vi\^) of the ith component. The 
unstable modes, which trigger CSI, should satisfy the condition 5J = {) with 
5J\ = —SJ2 due to the law of momentum conservation. Since the CSI is 
the dynamic instability triggered by unstable modes with an imaginary part 
Im w 7^ [45], the amplification of the unstable modes will exponentially 
enhance the momentum exchange. 

Figure 1.4 depicts the characteristic nonlinear dynamics of the CSI in an 
uniform two-component BEG obtained by numerically solving the GP equa- 
tion (1.7). In the early stage of the dynamics, amplification of the unstable 
modes creates disk-shaped low-density regions that are orientated in the x 
direction [Fig. 1.4(b)]. The lowest density inside the disk region reaches zero, 
creating a local dark soliton. The soliton in the ith component transforms 
into a vortex ring via snake instability [49] with a momentum antiparallel to 
the initial velocity Vi [Fig. 1.4(c)[. The vortex ring distribution can be deter- 
mined by the characteristic wavenumber of the unstable mode ~ toVr/Ti 
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(a)('=0 (bV'=12.2 (c)r'=12.6 (d) /'= 13.8 (eW'=26.0 (f)/'=34.1 




Fig. 1.4. Characteristic nonlinear dynamics of CSI of two-component BECs. The 
parameters are set to mi = m2 — m, gn = g22 = g, and g\2 ~ 0.9g, similar to those 
in experiments. In the numerical simulation, the initial state is prepared by adding 
small random noise to the stationary wave function •Z'io with nio = n2o ~ no and 
Vi — — V2. The panels shows the time development of low-density isosurfaces with 
ni(r) — 0.05no for Vi — ^Vrx. Because of the symmetric parameter setting, the 
second component behaves in a similar manner to the first component (not shown) 
[Takeuchi et al: Phys. Rev. Lett. 105, 205301 (2010), reproduced with permission. 
Copyright 2010 by the American Physical Society]. 

and k± ~ ■mVu/2h for a large relative velocity Vr > V+. The length scale 
~ k/Vr vifith K — 2Trh/m then characterizes both the radii of the vortex rings 
and the intervals between the rings along Vr immediately after vortex nucle- 
ation. Thus, the vortex line density ly after the instability is roughly estimated 
to be ly ~ V^/k"^, which can be controlled by varying the relative velocity Vr. 

Momentum exchange accelerates after vortex ring nucleation. The vortex 
motion then dominates the exchange. Since the momentum carried by a vortex 
ring increases with increasing radius, the radii of the nucleated vortex rings 
increase with time for momentum exchange. This dynamics resembles that 
of quantized vortices under thermal counterflow of liquid helium [50] , where 
the vortices are dragged by the mutual friction between the superfluid and 
normal fluid components. When the vortex rings become large, the interaction 
between the vortex rings deforms the rings and vortex reconnections occur, 
which depresses the momentum exchange [Fig. 1.4(d)]. These effects make the 
vortex dynamics very complicated, leading to binary QT in which the vortices 
of both components are tangled with each other [Fig. 1.4(e)]. The momentum 
exchange almost terminates and each component has an average momentum 
of nearly zero. 

The relative motion of two-component BECs and the resulting CSI can 
be experimentally realized by employing the Zeeman shift of atomic hyperfine 
states. Hamner et al. [48] initially prepared overlapping two-component BECs 
of *^Rb in the hyperfine states \F,mF) = |1,1) and |2,2), which satisfy the 
miscible condition 511322 > 9x2- When a magnetic field gradient was applied 
along the longer axis of the trap, the gradient generated forces in opposite 
directions for the two components due to the Zeeman shifts. Counter superflow 
then occurred and its relative velocity was controlled by the magnetic field 
gradient [48,51]. 
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Interface instability 

Next, we consider the interface instability of phase-separated two-component 
BECs. The interaction parameters satisfy the immiscible condition 511522 < 
§12 ■ We assume that components 1 and 2 arc phase separated at the interface 
near the y = plane, which is sustained by the external potential Ui{y). The 
density distributions ni{y) are also assumed to depend only on y and ni = 
for y > r] and n2 = for y < -q, where y = r]{x, z, t) is the interface position. 
Here, we neglect the interface thickness for simplicity. The Lagrangian of Eq. 
(1.4) can then be rewritten as 

L = Jdxdz(^J ' dyPi+ J dyP^ - aS, (1.13) 

where a is the interface tension coefficient a [52,53], which originates from the 
excess energy at the interface, and S = J dxdz[l + (drj/dx)^ + {d'qjdzf^^'^ is 
the interface area. Taking the functional derivative of the action J dtL with 
respect to 77(0;, 2, t) and setting it to zero, we obtain 

Pi (2/ = ^) - P^V = ^) + a + 0) = 0, (1.14) 

which corresponds to the Bernoulli equation in hydrodynamics. 

We consider a stationary state in which the ith component flows with a 
velocity Yi — ViX as tf'io = \/ni(j/)e*(™*^''^"^''"' which is similar to Eq. 
(1.8) except that the density depends on y. Substituting this into Eq. (1.14) 
with 7] = gives the equilibrium condition for the pressure, ^urii (0)^/2 = 

322^2(0)72. 

To analytically derive the dispersion relation for the interface wave, we 
assume that the system is approximately incompressible. We consider the 
small phase fluctuation iZ'i = iZ'joe"'^' and the interface mode r] = 67] as 

56i = Aie-^-'^^*^y cos{kx - ut), (1.15) 
6r] = asm{kx — ojt), (1-16) 

where Ai and a are infinitesimal parameters. From the kinematic boundary 
condition, the interface velocity in the y direction (d/dt + Vid/dx)rf must be 
equal to Ti/ {imini)'^* d'Pi/ dy\y=ri, giving 

-i-iy—Aike-^-^y'''^ = (Vik - uj)a. (1.17) 

mi 

Substituting Eqs. (1.15)-(1.17) into Eq. (1.14) and neglecting second and 
higher orders of Ai and a, we obtain 

^{uj - Vikf - /in,i = -y - V2kf - hn,2 + ak\ (1.18) 
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where n^i = ni(r? - 0), 71^2 = "2(77 + 0), pj = miTisi, and /j = dUi/dy\ 
Equation (1.18) gives the dispersion relation, 



where Vq = {piVi + /02V2)/(pi + P2) is the center-of-mass velocity and F = 
nsifi — ns2f2 is the force due to the gradient of the external potential. We note 
that Eq. (1.19) has the same form as the dispersion relation for an interface 
wave in classical incompressible and inviscid fluids. In fluid dynamics, the 
gradient of the potential is equivalent to gravity. 

Rayleigh-Taylor instability 

First, we consider the case Vi = V2 = 0- For F < 0, the system is always 
dynamically unstable in the wavenumber range < fc < y^\F\/a, which is 
known as a Rayleigh-Taylor instability [54, 55] . This situation corresponds to 
a layer of a lighter fluid under a heavier fluid layer in a classical fluid, where 
the translation symmetry of the interface is spontaneously broken. 

Sasaki et al. [54] proposed a system of two immiscible BECs with different 
hyperfine spins (e.g., \F,niF} = |1,— 1) and |1,1) of ^'^Rb atoms) placed in 
an external magnetic field gradient B' = dB/dz. Such condensates experience 
the potentials +iibB' z/2 and where jiB is the Bohr magneton. 

The force generated by this potential can realize F < so that the two 
condensates are pushed in opposite direction. Numerical simulations of the GP 
equations reveal that this gradient modulates the interface so that it grows in 
a mushroom pattern. Vortex rings then nucleate due to atoms near the center 
flowing upward and atoms at the periphery of the cap of the mushroom shape 
flowing downward. Gautam and Angom [55] considered a system of a ®^Rb- 
^''Rb BEG mixture and the Rayleigh-Taylor instability caused by tuning the 
interspecies interaction through a Feshbach resonance. The signature of the 
instability should appear in the damping behavior of the collective shape 
oscillation. 

Richtmyer-Meshkov instability 

The Richtmyer-Meshkov instability occurs when an interface between fluids 
with different densities is impulsively accelerated (e.g., by the passage of a 
shock wave). For atomic BEGs, this instability can be caused by a magnetic 
field gradient pulse B'{t) oc S{t) [56]. The nonlinear stage of this evolution 
is qualitatively similar to that of the Rayleigh-Taylor instability. However, 
the instability dynamics differs considerably from that of classical fluids. The 
main difference originates from the quantum surface tension and capillary 
waves, which suppress perturbation growth and droplet detachment from an 
elongated perturbation finger. The instability for more general time-dependent 
forces has been discussed in Ref. [57] 




(1.19) 
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Kelvin-Helmholtz instability 

Finally, we consider a system with the shear flow Vi 7^ V2 7^ (Vr ^ 0). The 
dispersion relation Eq. (1.19) implies that, for > 2\/Fa{pi + P2)/ P1P2 = 
V[|jj, the imaginary part Im(u;) becomes nonzero and the shear-flow states are 
dynamically unstable against excitation of the interface modes with k- < k < 
k+ , as in classical Kelvin-Helmholtz instability, where k± — kozL \/k'^ — F/a 
with fco = pip2V'^/2a{pi + p-i). When F = 0, ^kh vanishes and the system is 
always dynamically unstable for > 0. In addition to dynamic instability, 
thermodynamic instability can occur due to dissipation when < [58] . Here, 
we restrict ourselves to dynamic instability in nondissipative systems. 




Fig. 1.5. Nonlinear dynamics of Kelvin-Helmholtz instability in phase-separated 
two-component BECs with a relative velocity Vn = Vi — V2 > Vkh. The numerical 
simulation of the GP equation was performed with a period L = 64^ along the a;-axis 
and the Neumann boundary condition with L — 32^ in the j/-direction, where ^ is 
the healing length. The unit of time is ^/c, where c is the sound velocity. The initial 
state (a) is with a small random seed to trigger the instability. (Upper panels) 
The height and color show the vorticity cjca = (V x Vg)z and the density difference 
ni — n2 between the two condensates, respectively. (Lower panels) Two-dimensional 
plots of ni -712. [Takeuchi et al: Phys. Rev. B 81, 094517 (2010), reproduced with 
permission. Copyright 2010 by the American Physical Society]. 

Figure 1.5 demonstrates the Kelvin-Helmholtz instability for Vr > Vkh 
[58] . In the linear stage of the instability, the sine wave corresponding to the 
most unstable mode with the maximum imaginary part max^ {Im [w(/c)]} is 
predominantly amplified. As the amplitude increases, the sine wave is dis- 
torted by nonlinearity [Fig. 1.5(b)], and deforms into a sawtooth wave [Fig. 
1.5(c)]. The vorticity u>cs increases on the edges of the sawtooth waves and cre- 
ates singular peaks [Fig. 1.5(d)]. Subsequently, each singular peak is released 
into each bulk, becoming a singly quantized vortex [Fig. 1.5(e)]. The release 
of vortices reduces the vorticity of the vortex sheet and therefore reduces the 
relative velocity across the interface. The released vortices drift along the in- 
terface and the system never recovers its initial flat interface. These nonlinear 
dynamics differ considerably from those in classical KHI, where the interface 
wave grows into roll-up patterns. 

The above discussion is valid when the interface thickness ^ 912/9 — 1 
(for the case 511 — 522 = 9) is much smaller than the wavelength of the unsta- 



16 Makoto Tsubota^ and Kenichi Kasamatsu^ 

ble interface mode, typically given by ~ U/mVn. In the opposite case, the CSI 
becomes the dominant instability of the flowing state. The crossover relative 
velocity between the two instabilities is evaluated as ~ hy/ 912/ g — 1/m^ 
[59]. 



I. 4 Conclusions 

Since the realization of BEC in a dilute atomic gas in 1995, most studies on 
its QHD have been limited to vortex lattices under rotation or motion of a few 
vortices. Howciver, as research on superfluid helium has shown, there are many 
other interesting problems on QHD in atomic BECs, some of which have been 
discussed in this paper. Important topics are quantum hydrodynamic instabil- 
ity and QT beyond this instability. Unlike classical fluid dynamics, most QHD 
phenomena can be reduced to the motion of quantized vortices. For example, 
for QT in atomic BECs, the cascade process of quantized vortices, which 
transfer energy from large to small scales, can be visualized. The observa- 
tion of Kolmogorov spectra could confirm the cascade process in wavenumber 
space. The observation of QT in atomic BECs enables us to combine cascade 
processes in real and wavenumber spaces. Such investigations are almost im- 
possible in superfluid helium and in classical turbulence. We anticipate that 
this research field will develop rapidly in the near future. 
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